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Abstract 

Recent observations from type la Supernovae and from cosmic microwave background (CMB) anisotropics 
have revealed that most of the matter of the Universe interacts in a repulsive manner, composing the so- 
called dark energy constituent of the Universe. The analysis of cosmic gravitational waves (GW) represents, 
besides the CMB temperature and polarization anisotropics, an additional approach in the determination of 
parameters that may constrain the dark energy models and their consistence. In recent work, a generalized 
Chaplygin gas model was considered in a flat universe and the corresponding spectrum of gravitational waves 
was obtained. The present work adds a massless gas component to that model and the new spectrum is 
compared to the previous one. The Chaplygin gas is also used to simulate a A-CDM model by means of 
a particular combination of parameters so that the Chaplygin gas and the A-CDM models can be easily 
distinguished in the theoretical scenarios here established. The lack of direct observational data is partialy 
solved when the signature of the GW on the CMB spectra is determined. 
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1 Introduction 

A large number of dark energy candidate models ha ve been proposed sinc e that the SNela |Riess et all Il998l 
IPermultter et all Il998l iTonrv et all [2003] and CMB |Tegmark et all l2004| experiments revealed the accelerat- 
ing expansion of the universe. In order to explain the observational data all these models consider an equation 
of state with negative pressure. The cosmological constant model s consi der a simple equation of state, p = — p, 
which however results a huge discrepancy wit h the data |Sahnil . |2004j . Another class of models considers a 
scalar field, called quintessence |Caldwell l2004j . a third class assumes a perf ect fluid with a neg ative pressure 
which is proportional to the inverse of the energy density, the Chaplygin gas |Gorini et all l2004j . and finally a 
forth class, considers an X- flu id with an e quation of state p = —uop, where a; is a positive constant, known as 
phantom energy when uj > 1 Salmi l2004j . 

The distinction among these models must be done by means of combined observations which up to the 
present have been unable to define what description is more appropriate. Gravitational waves represent one of 
the potential tools that may in near future offer an additional way to constrain the cosmological parameters 
and distinguish among the models, since gravitons decouple in a very early time in the universe history. In this 
paper, we study the spectrum of gravitational waves due to the X-fluid model and compare the results with a 
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previous work |Fabris et all l2004| , where the cosmological constant and the generalized Chaplygin gas models 
were taken into account. In section ^3 these models are described; in J3the GW equations are presented. The 
spectra analysis is left to which is followed by the conclusions of this work. 



2 Outline of the model 

We consider a flat, homogeneous and isotropic universe described by the Friedman- Robertson- Walker metric, 
which makes the Einstein's equations to assume the form 



(Pm + Pde), 



(1) 
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where a is the scale factor of the universe, while p m and pde are the pressureless fluid and the dark energy 
densities, respectively. The pressures p m and pd e of the fluids are related with their densities by the equations 
of state p m = 0, and p x = cup x , with uo < (in case of an X-ffuid) or p c = —A/p^, with A, a > (in case of 
Chaplygin gas), respectively. We take the scale factor today as the unity, ao = 1 (The subscripts 0, according 
to the current notation, indicate the present values of the quantities.) and rewrite the dark energy density for 
both cases as: 
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With equations J3J and 
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we write the Einstein's equations in the form 
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for the Chaplygin gas. The Hubble constant Ho is defined as Ho = do /a and, once we are restricted to a flat 
universe, the fractions of pressureless matter and dark energy gas today, ^ mo and flde , respectively, obey the 



relation ^ mo + Q deo = 1. 

With these two last equations we are able to write the GW amplitude differential equation as a function of 
the observable parameters Hq, a, ^ mo and O^ eo , and of the dark energy fluid parameters: A, a or uo. It is also 
important to remark that 

(i) if A = 0, then the Chaplygin gas behaves like the pressureless fluid and the situation is the same as if we 
had set ll mo = 1; 

(ii) on the other hand, for A = 1 it behaves like the Cosmological Constant and therefore, we can simulate 
the A-CDM scenario; 

(iii) the X- fluid is equal to the Cosmological Constant for uo = — 1. 

Among the many possible cases produced by the combinations of the parameters we have chosen a few 
important ones which are indicated in table 1. These choices will be sufficient for the aim of the present 
analysis, where we compare the spectra due to each model. 
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Table 1: Models chosen for the analysis. 



3 GW equations 

Cosmological gravitational waves are obtained through little perturbations hij on the metric. Hence, the tensor 
Q^j \ related to the unperturbed metric, is replaced by g^ = gff +hij and the resulting expression is | Weinberg . 

h--h+(^-2-)h = 0, (9) 
a \a z a J 

where k is the wave number times the velocity of light (k = 27rc/A), the dots indicate time derivatives and 
h(t) is defined as: hij(x,t) = h(t)Qij, where Qij are the eigenmodes of the Lagrangian operator, such that 
Qu = dkQki = 0. 

Proceeding with a variable transformation from time to the scale factor a, and representing the derivatives 
with respect to a by primes, equation ® assumes the form 

h'> + (±-±)h' + ±(£-2*)h = Q (10) 
\a z a J a z \a z a J 

and, by the application of the background equations J7J) and 0, or JSJ) and ©, which concern to the dark 
energy fluid, into (fTU|) . one can easily express h in terms of the parameters of the model. 



3.1 Generalized Chaplygin gas 



Let us perform the last operation mentioned above and, in order to find h as a function of the redshift z, let us 
(i) use the known relations l + z=^,ao = l; (h) perform a second change of variables (from a to z)\ and (hi) 
take back the dots to indicate, from now on, the new integration variable. These operations result in 
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where k is redefined to absorb the Hubble constant, i.e., k = 2ttc/Ho\ and, therefore, h is a dimensionless 
quantitiy. 
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3.2 X-fluid 

To obtain the equivalent equations for the X-fluid model, we use the same procedure described above. The 
result may be written in the same form as in l|1 1|> . with /i and f 2 redefined as: 



/i = n mo + n xo (i + w )(i + z) 
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f 2 = n mo (i + z) 3 + n X0 (i + z) 
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One can easily notice that, if uj = — 1 and A = 1, both cases converge to the same equation and then reproduce 
the A-CDM scenario. 



4 GW spectra analysis 

The power spectrum of gravitational waves, defined as 



dhiv 



= \h Q {v)\v h ' 2 



(16) 



[where ho{v) = h(0) and v = Hok/27t] : is generally obtained from the solutions of {TT} which are found by means 
of a numerical algorithm specifically created for this problem. In the particular case of the A-CDM model, an 
analytical result is also possible and this fact is used to verify the accuracy of the calculation. 

For each of the cases of interest, we have assigned some common parameters, namely the initial conditions 
h(zi) = ^" 3 10" 5 , h(zi) = ^" 2 10" 5 , Zi = 4000; the range of frequencies 10" 18 i^2: <v< 10~ 15 Hz, which is a 
small part of the range of obser vational interest for GW that goes up to up to 10 10 i^2:; and the normalization 
constant imposed by the CMB |Riazuelo fc Uzanl . Ec 
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The resulting spectra are presented in figures (1-11) and we can see that they all have a similar behavior: 
a strong oscillatory shape; a peak at v ~ (due to the initial condition, which is proportional to /c -3 ); after a 
very rapid decrease, the spectra increase slowly. Besides these features, it is important to remark that: 
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Figure 1. GW spectrum for Q CQ = and flm = 1. 
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Figure 2. GW spectrum for A = 1, Q CQ = 0.96 and Qm Q = Figure 3. GW spectrum for A = 1, Q CQ = 0.7 and Qm Q 

0.04. 0.3. 



(i) The plot (Fig. 1) is for a pressureless pure model and we can see that the spectrum corresponding to 
a Cosmological Constant dominated universe (Fig. 2) grows faster with the frequency reaching a value three 
times greater. Still referring to a A-CDM model, figure (Fig. 3) shows a behavior similar to the pressureless 
fluid, but the growing amplitude is just a little larger. 
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Figure 4. GW spectrum for A = 0.5, a = 1, Q CQ = 0.96 and Figure 5. GW spectrum for A = 0.5, a = 1, Q CQ = 0.7 and 



Qm = 0.04. 



^mn — 0.3. 
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Figure 6. GW spectrum for A = 0.5, a = 0.5, Q CQ = 0.96 Figure 7. GW spectrum for A = 0.5, a = 0.5, Q CQ = 0.7 and 

and Q mo = 0.04. ft mo = 0.3. 
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Figure 8. GW spectrum for A = 0.5, a = 0, Q CQ = 0.96 and Figure 9. GW spectrum for A = 0.5, a = 0, ^ CQ = 0.7 and 
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- 0.3. 



(ii) The graphics (4-9), however, are not so easily distinguishable. They seem to have no dependence on 
the Q Co parameter and this may be interpreted as a consequence of the Chaplygin gas property of interpolation 
between dark-matter and dark-energy. 

(hi) Figures (4) and (5) refer to the Chaplygin gas in its traditional form - with the equation of state 
p = —A/p, since a = 1 - while (6-9) correspond to the generalized Chapligyn gas and the six plots together 
show that the spectrum is not affected by a. 
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Figure 10. GW spectrum for uj = —0.5, Q x 

^mn — 0.3. 
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(iv) The X- fluid is represented in plots (10) and (11). The later is the case of phantom energy, with uj = — 10, 
while the former is an intermediate situation between the pressureless and the A cases. We notice that they 
also do not differ from each other. 



5 Conclusions 

In this work we have compared three important dark energy models in the context of the gravitational wave 
spectra that they are able to induce and we observe that the models are strongly degenerated in the range 
of frequencies studied - from 10 -18 up to 10~ 15 Hz. This degeneracy is in part expected, since the models 
must converge to each other when some particular combinations of parameters are considered. In addition, 
for the Chaplygin gas, t he negligible dependence on the density parameter ^ Co is consistent with the idea of a 
unified dark component |Fabris et all 1200*3. iGorini et al.l . l2004j . Since the gas interpolates both dark energy and 
matter, the suppression of the dark matter density parameter (see, plots (4-9)) should not affect the spectrum, 
as observed. However, it is interesting how it seems to be insensitive to a as well. This result indicates that this 
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model is more sensitive to A than to a. The X-fluid component contribution seems to be completely degenerated 
and must be better studied in the near future, specially at larger frequency ranges. 

Unfortunately, no direct observational data from cosmological gravitational waves is available up to now, 
but still it is important to investigate the relation between these gravitational waves produced or amplified 
by the presence of the dark energy co mponent and the CMB temperature and polarization anisotropics which 
are already observationally determined [Bennett et all l2003j . This correlation between CMB and gravitational 
waves will be properly studied once we have a larger range of frequencies in the GW spectra. 
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